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Abstract 



The decay rate of the r lepton into hadrons of invariant mass smaller 
than Q 3> Aqcd can be calculated in QCD using the OPE. Using experi- 
^ ■ mental data on the hadronic mass distribution, the running coupling con- 

stant a s (Q 2 ) is extracted in the range 0.85 GeV < Q < m T , where its value 
changes by about a factor 2. At Q = m T , the result is a s (m 2 ) = 0.33±0.03, 
corresponding to ot s {m 2 z ) = 0.119 ± 0.004. The running of the coupling 
constant is in excellent agreement with the QCD prediction based on the 
three-loop /3-function. 
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1 Introduction 



One of the most accurate methods to determine a s in the low-energy region is 
provided by the measurement of R T , the r decay rate into hadrons normalized to 
the leptonic decay rate: 

r(r — > v T + hadrons) , . 

i? r can be calculated in QCD using the Operator Product Expansion (OPE) 
[HI 01. The result is: 



R T = N c jl + 5 pcrt [a s (m 2 T )} + 5 power } 



nJi + + 5.202 f^SV + 26.37 ' 



I 7T \ 7T / \ 7T 

_8|K/4 + 32^<^-2M + ...}. (2) 

The non-perturbative power corrections in this expression are proportional to the 
strange-quark mass, the quark condensate, and higher- dimensional condensates. 
Because all contributions of dimension less than six vanish in the chiral limit, 
the power corrections are numerically small; using standard values of the QCD 
parameters, one finds £ P ower = —(1-4 ± 0.5)%. This, together with the fact that 
the perturbation series is known to third order, make R T a good observable to 
measure a s . 

Experimentally, R T is obtained from the relation R T = 1/B e — 1.97256, where 
B e is the leptonic branching ratio. Direct measurements give B e = (17.80 ± 
0.06)% 0, whereas using the r lifetime, r T = (291.3 ± 1.6) fs g], we obtain 
B e = t t /t^ (m T /m^) 5 = (17.84 ± 0.10)%. Averaging the two results gives R T = 
3.642 ±0.010, and taking into account small electroweak radiative corrections not 
displayed in @ we obtain 

5 pcrt [a s (m 2 T )} = 0.205 ± 0.003 cxp ± 0.005 th . (3) 

The dominant theoretical uncertainty in extracting a s from this measurement 
comes from the truncation of perturbation theory || , which induces an error of 
order a 4 s . This uncertainty can be estimated by considering some approximate 
resummations of the perturbation series (starting at order a\) and comparing 
them to the fixed-order calculation. The resummation procedure of Le Diberder 
and Pich || resums certain "large- 7r 2 " terms to all orders in perturbation theory. 
Recently, another class of terms, the so-called renormalon chains (J7|, have been 
investigated. These are the terms of order /3o _1 a™ in the perturbation series for 
^pcrt) where (3q is the first coefficient of the QCD /5-function. The resummation 
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Fig. 1. Different perturbative approximations for the quantity <5 per t: exact 
order-a^ result in the MS scheme (solid line), resummation of renormalon 
chains (dashed line), resummation of Le Diberder and Pich (dash-dotted 
line). The experimental result (||) is shown as a band. 



of such terms in the case of R T has been discussed in Refs. [|[ §. In Fig. [T], we 
show the corresponding theoretical predictions for <5 pert as a function of a s (m 2 ). 
We conclude that Sa s (m 2 ) ~ ±0.03 is a reasonable estimate of the truncation 
error. This leads to 



a Am 



I) = 0.33 ±0.03, a s (m%) = 0.119 ±0.004. (4) 



For the sake of completeness, we have translated our result into a value of a s at 
the mass of the Z boson. 

The analysis just described provides one of the best determination of the 
QCD coupling constant in the low-energy region. The result (§) is included 
in Fig. ^ which shows a collection of measurements of a s performed at different 
energy scales fllQf . Besides r decays, low-energy (Q ~ 1.6-10 GeV) measurements 



come from deep-inelastic scattering and T spectroscopy and decays. At higher 
energies (Q ~ 30-130 GeV), the most reliable determinations of a s come from 
measurements of the total cross section, jet rates and event shapes in e + e~, pp 
and ep collisions. Taken all together, these measurements provide clear evidence 
for the "running" of the effective coupling constant a s (Q 2 ), which in QCD is 
predicted to decrease with the momentum transfer. This property of "asymptotic 
freedom" |ll| is one of the key predictions of QCD. Formally, it is expressed by 



the fact that the /^-function is positive, where 
^1 = -a s (Q 2 ) /?MQ 2 )] 



2 



(5) 



and Pq = 9, (3\ = 64 and Pi — 3863/6 are the first three expansion coefficients 
of the /3-function, evaluated for rif = 3 light quark flavours. (The value of Pi is 
specific to the MS renormalization scheme.) 
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Fig. 2. Compilation of a s measurements. The curves correspond to the QCD 
prediction for the running coupling constant for a s (m 2 z ) = 0.116 ± 0.005. 



A test of the running of a s by combining measurements performed in many 
experiments operating at different energy scales has the disadvantage of involving 
different experimental systematic errors, as well as different levels of sophistica- 
tion of the theoretical calculations. Therefore, it is an appealing idea to measure 
the scale dependence of the coupling constant in a single experiment. This can be 
done in high-energy experiments at pp and ep colliders, where a large range of Q 
values can be probed simultaneously ||12|| . However, so far the precision obtained 
in these measurements is rather low. In this talk, we propose a high-precision 
test of the running of a s in the low-energy region (0.85 GeV < Q < m T ), using 
data obtained in a single experiment |13| . The value of a s changes by about a 
factor 2 in this energy range, which is equivalent to the variation between 5 and 
100 GeV. 
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2 Extraction of cx s (Q 2 ) in r decays 



We shall consider the r decay rate into hadrons of invariant mass squared smaller 
than so, normalized to the leptonic decay rate: 

r(r -> i/ T + hadrons; s ha d < s ) /" , di? T (s) 

so = ^ = / ds — , 6 

1 r — > v T e v e ) j ds 

o 

where di? r /ds is the inclusive hadronic spectrum, which has been measured by 
the CLEO and ALEPH Collaborations 0, HI. To obtain dR T /ds, we have 



multiplied the normalized distributions by R T . The result is shown in the upper 
portion of Fig. [| Not shown in the figure is the contribution from r — > h~v T with 
h~ = 7r~ or K~ , which has a branching ratio of (11.77 ± 0.14)% |J. Integrating 
these spectra over s and combining the results weighted by their statistical errors, 
we obtain the distribution R t (sq) shown in the lower portion of the figure. Sys- 
tematic errors have been estimated by taking the difference between the CLEO 
and ALEPH data, and added in quadrature with the statistical errors. Since the 
errors are strongly correlated, the result is presented as a band. 

Using the analyticity properties of QCD spectral functions, the quantity 
R T (s ) can be represented as a contour integral along a circle of radius |s| = s 
in the complex s-plane (for simplicity, we quote the result in the chiral limit): 

Rr{* ) = ± I -w(*L,^)d(s). (7) 
2m J s \mi mi J 

ii ' ' 

\s\=so 

Here 

w(x, y) = 2{x -y)- 2(x 3 - y 3 ) + (x 4 - y 4 ) (8) 

is the phase-space function, and D(s) is a current-current correlation function, 
which contains all QCD dynamics. The representation (^) shows that so is the 
only scale at which QCD dynamics is probed; the r-lepton mass appears only 
in the phase space. Provided that s ^> Aq CD , the correlation function D(s) 
is needed at large momentum transfer only, and the OPE can be employed to 
calculate -R r (s ) as a function of a s (s ) and x = s /m^: 

R T {s ) = N c {r pcrt [a s (so),£o] + ' power (so)}- (9) 
The perturbative contribution is given by (oq = a s (so)/7r): 



pe rt[a s (so),x ] = (2x - 2x1 + x o) 1 + a o + lM0a 2 - 10. 28^ + (K 4 - 156. 0)a 
+ (2x - \xl + \xt) (2.25ag + 11.38ag - 46.24a^) 
+ (2x - lx 3 + ±xl) (10.125ag + 94.81a 4 ) 
+ (2x - A x 3 + ^^4) 6g 34a 4 + 0(a 5 } ) (1Q) 
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Fig. 3. Upper plot: The hadronic mass distribution di? T /ds in r decays. 
Shown are the statistical errors after correcting for detector effects. CLEO 
data from electron-tagged and muon-tagged events have been combined. 
The ALEPH data are preliminary. Lower plot: The integrated spectrum 
R T (so). The experimental result, including statistical and systematic errors, 
is presented as a band. The curves will be explained in Sect. 0. 



where K± is the five-loop coefficient in the Adler function, which is currently not 
known exactly. In our analysis, we use the estimate K$ ~ 27.5 [TH] obtained 



using the methods of Ref. [17]. As in the case of R T , the truncation of the 



perturbation series will turn out to be the main theoretical uncertainty in our 
analysis. We estimate the importance of the unknown higher-order contributions 
(of order Oq and higher) by resumming the renormalon-chain contributions to all 
orders in perturbation theory, using the results of Ref. 0. We shall compare 
fixed-order perturbation theory with this resummation and take the difference 
as an estimate of the perturbative uncertainty. This estimate of the truncation 
error is more conservative than that obtained by dropping the last term in the 
series in ([Top. 
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Fig. 4. Values of a s (so) extracted from the data on R t (sq). The dark band 
represents the experimental errors, the light one the sum of the experimental 
and theoretical errors. The errors are strongly correlated. The dashed line 
shows the three-loop QCD prediction for the running coupling constant. 



The power corrections in ([|) are given by 
r p0 w Cr (so) = -6 \V US \ 2 (l + x -xl + lag) 4 + 32vr 2 _ 2 M + .... (1 1) 

Note that, as a simple consequence of the representation (|7p, no inverse powers 



of So appear in this expression This is true as long as the coefficients of the 
power corrections to the correlation function D(s) do not contain logarithms of s. 
As a result, the OPE converges well down to low scales sq- For instance, we find 
fpower = —(1-4 ± 0.5)% at s = ml, and r powcr = -(1.5 ± 0.5)% at s — 1 GeV 2 . 
The break-down of the OPE (see Sect. [3] below) will thus not be driven by a 
blow-up of the series of power corrections. Another important feature of (|TT| ) is 
that the terms involving the vacuum condensates are independent of So- Hence, 
the uncertainties in the values of the condensates do not affect the so dependence 
of R T (so), which will be used to study the running of a s (so). 

From the measurement of the quantity R t (sq) shown in Fig. |3], we extract 
a s (s ) as a function of s by fitting to the data the theoretical prediction ob- 
tained using fixed-order perturbation theory. The result, including experimental 
errors only, is represented by the dark band in Fig. |j. Theoretical uncertainties 
arise from the truncation of the perturbation series and from the uncertainty in 
the values of the nonperturbative parameters. They affect the overall scale of 
the a s values (by about 8-10%), but have very little effect on the evolution of 
the coupling constant. The sum of the experimental and theoretical errors is 
represented by the light band. The dashed curve shows the QCD predictions for 
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a s (so) obtained at three- loop order, normalized to the central value of the data 
at So = vn? T . The observed scale dependence of the running coupling constant is 
in excellent agreement with the QCD prediction. 




Fig. 5. Values of a s (so) extracted from the data on R T (so) using fixed-order 
(FOPT) and resummed perturbative theory (RPT). The dashed lines show 
the QCD prediction obtained using the three-loop /3-function. The dash- 
dotted lines refer to the one-loop /3-function. 

As mentioned above, the main theoretical uncertainty comes from the trunca- 
tion of the perturbation series. This is illustrated in Fig. ||, where the evolution 
of a s (so) as a function of sq is shown separately for fixed-order and resummed 
perturbation theory. The curves show the QCD predictions for the running cou- 
pling constant obtained at one- and three-loop order, normalized to the data at 
So = vn? T . It is seen that higher-order corrections effectively renormalize the over- 
all scale of the a s values (i.e. the Aqcd parameter). For instance, the value of a s 
at So = vr? T changes from 0.33 (fixed-order) to 0.31 (resummed). The difference 
between the two results for a s (so) has been used to estimate the truncation error. 

Fig. |6| shows our result combined with the other measurements of a s collected 
in Fig. [| We have replaced the data point at Q — m T by the band shown in 
Fig. f|, which extends to much lower values of Q. This figure demonstrates nicely 
the main features of our approach: it extends the range of a s values accessible 
to experiments, thus allowing a measurement of the strong coupling constant at 
scales lower than the lowest ones attainable before. Moreover, it provides a test of 
the QCD evolution of a s with higher precision than all other single measurements 
of the running to date. 
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Fig. 6. Compilation of a s measurements including our result obtained from 
the analysis of hadronic r decays. The curves correspond to the QCD pre- 
diction for the running coupling constant for a s (m 2 z ) = 0.121, 0.116, 0.111 
(top to bottom). 



3 Break-down of the OPE and quark— hadron 
duality 

An important question which we have to address is to determine the lowest value 
of So for which our analysis can be trusted. In other words, at which point do we 
expect the OPE to break down? To answer this question, it is important to realize 
that we are applying the OPE in the physical region (i.e. the region of time-like 
momenta), where QCD cannot be used to calculate correlation functions such as 
D(s). The reason why we trust the calculation of R t (sq) is that to perform the 
contour integral in (0) requires knowledge of the correlation function for large 
(complex) momenta only. Moreover, the integrand vanishes for s = sq, where the 
contour touches the branch cut of -D(s); hence, the main contributions come from 
regions far away from the singularities, where the OPE can be applied. Another 
way to say this is that in the calculation of R T {so) we assume quark-hadron 
duality, which is the hypothesis that QCD can be employed to calculate physical 
decay rates if they are "smeared" over a sufficiently wide energy interval [|TE|]. 
In the present case, this smearing is provided by the integration over the range 
< s < so in (|5|). The question of how accurate the duality assumption is and for 
what values of sq it applies is, however, a phenomenological one. Despite of some 
interesting new ideas [plU , it cannot be answered yet from theoretical grounds. 
To test the assumption of duality, we compare the data for the quantity 
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Rt(so) with the theoretical predictions obtained from the OPE, using both fixed- 
order and resummed perturbation theory. The results are shown by the two 
curves in the lower portion of Fig. |3|. In obtaining these curves, we have adjusted 
the value of a s (m 2 ) so as to fit the data at so = m 2 . The value of a s (so) 
is then obtained from the solution of the renormalization-group equation (^). 
Theoretical uncertainties have little influence on the Sq dependence of R t (sq). For 
the perturbative part of the calculation, this is apparent from the good agreement 
of the two theoretical curves in Fig. |3|, which refer to values of a s (m 2 ) that differ 
by 9%. Hence, the so dependence of R T {so) is predicted essentially without any 
free parameters, and the comparison of the data with the theoretical predictions 
provides a direct test of quark-hadron duality. 

We find excellent agreement over the range 0.7 GeV 2 < s < m 2 , indicating 
that in r decays duality holds as soon as the integral over the hadronic mass 
distribution includes the p resonance peak. This justifies a posteriori our choice of 
the energy interval in the previous section. It is remarkable that, once sq exceeds 
the value of 0.7 GeV 2 , the onset of duality happens almost instantaneously. Since 
the p meson is such a prominent resonance, this is the best possible scenario that 
could be expected. The small oscillation of the experimental band around the 
theoretical curve, which could be due to some deviations from duality in the 
region of the a\ resonance, are not significant given the precision of the data. 
Even if such oscillations will be confirmed in further analyses based on more 
precise data, they will clearly not put a severe limitation on the applicability of 
our method. 



4 Measurement of the /3-function 

To quantify the agreement between the data and the QCD prediction for the 
running coupling constant exhibited in Figs. ^ and [|, we extract from the data 
the /3-function defined in (J^j and compare the result to the prediction of QCD 
perturbation theory. Introducing the variable x = a s (so)/47r, we have 

4?r da s {s ) (3{x) 2 

p + Pix + (5 2 x + ... . (12) 



a 2 (s ) din s x 

We approximate the derivative da s /dlnso by a ratio of differences, Aa s /Alnso, 
for a selected set of so values chosen such that the differences Aa s are large 
enough to be significant given the errors in the measurement. For a s (so) in (|I~2j) 
we take the central value of each interval. We use the following so values: 0.75, 
0.95, 1.35, 2.06, and 3.16 GeV 2 , corresponding to four intervals of increasing 
width A In So, but constant Aa s ~ 0.075. The results are shown in Fig. [7]. The 
circles are obtained using fixed-order perturbation theory, while the squares re- 
fer to resummed perturbation theory. As expected, the two methods give very 
similar results for the running of the coupling constant. The estimate of the 
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Fig. 7. Experimental determination of the /3-function. The circles are ob- 
tained using fixed-order perturbation theory, the squares refer to resummed 
perturbation theory. The curves show the QCD /3-function at one-loop 
(dash-dotted), two-loop (dashed) and three- loop (solid) order. 



errors includes the theoretical uncertainties, the error due to the choice of finite 
intervals in ac s , and the experimental errors, which in this case are the dominant 
ones. The curves in Fig. |7| show the QCD /3-function at one-, two- and three- loop 
order in perturbation theory. The data provide clear evidence for the running 
of the coupling constant. Moreover, they prefer a running that is stronger than 
predicted at one-loop order. Indeed, between the three curves, the best descrip- 
tion of the data is provided by the three-loop prediction. Performing a fit with 
the three-loop /3-function, where (3 = 9 and (3i = 64 are kept fixed but the 
three-loop coefficient (3 2 is treated as a parameter, we find P^/Pf 1 = 1-6 ± 0.7 
using fixed-order perturbation theory, and 1.8±0.8 using resummed perturbation 
theory. 

We believe that such an experimental determination of the /3-function beyond 
the leading order can at present be done only in r decays. (A high-precision 
measurement of R e + e -(s) in the region below the charmonium resonances would 
provide an alternative place for such a study.) At higher energies, the value of 
a s is too small to distinguish between the three curves in Fig. 0; measurements 
in the region Q ~ 100 GeV, for instance, correspond to values x ~ 0.01. 

5 Conclusions 

We have presented a method to measure the running coupling constant a s (Q 2 ) 
in the low-energy region 0.85 GeV < Q < m T , using r-decay data obtained in a 
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single experiment. At Q — m T , we obtain a s (m 2 ) = 0.33 ±0.03, corresponding to 
the rather precise value a s (m 2 z ) = 0.119 ± 0.004. Our method provides a test of 
the scale dependence of the coupling constant in a region where this effect is most 
pronounced. The theoretical analysis is based on the OPE and the assumption 
of quark-hadron duality. We have tested this assumption and find that it holds 
provided the r decay rate is integrated over an energy interval large enough 
to include the p resonance peak. Our analysis provides a test of QCD at scales 
comparable with the lowest ones attainable before (Q ~ 1.6 GeV in deep-inelastic 
scattering), and with higher precision than all other single measurements of the 
running to date. We have extracted for the first time the /3-function from data 
and find that it is in good agreement with the three-loop prediction of QCD. 
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